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Abstract 

The Crab Pulsar is a relatively young neutron star. The pulsar is the central star in the Crab Nebula, a remnant of 
the supernova SN 1054, which was observed on Earth in the year 1054. The Crab Pulsar has been extensively observed 
in the gamma-ray energy band by the Large Area Telescope (LAT), the main instrument onboard the Fermi Gamma-ray 
Q_|Space Telescope, during its first months of data taking. The LAT data have been used to reconstruct the fluxes and the 
energy spectra of the pulsed gamma-ray component and of the gamma-rays from the nebula. The results on the pulsed 
component are in good agreement with the previous measurement from EGRET, while the results on the nebula are 
consistent with the observations from Earth based telescopes. 
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The Fermi Gamma- Ray Space Telescope was launched 
■by NASA on June 11, 2008 from Cape Canaveral, and 
started taking data in sky survey mode and in nominal 
science configuration on August 4, 2008. The main 
instrument onboard Fermi is the Large Area Telescope 
.(LAT) [H , a pair conversion telescope that covers a photon 
energy range extending from 20 MeV beyond 300 GeV. 
Fermi is moving on a quasi-circular orbit at 565 km 
'altitude with an inclination of 25.6°, and observes the 
.whole sky every 2 orbits (corresponding to 3 hours). 

The Crab Pulsar and nebula are the remnants of the 
explosion of a supernova reported by Chinese astronomers 
'in 1054. The Crab Pulsar, with a spin-down power 
.E = 4.6 ■ 10^^ erg s~^, is one of the most energetic known 
pulsars. It is located at a distance of about 2 kpc Q from 
Earth and is spinning with a period of 33 ms j3| . 

Both the Pulsar and Nebula have been extensively 
observed in all the energy bands, from radio to high energy 
'7-rays. In this paper we report the results of a spectral 
analysis performed on both the Crab Pulsar and Nebula 
using the data collected by the Fermi LAT during its first 
8 months of operation. 

2. Radio timing analysis 

The timing solution for the Crab Pulsar is built 
using the data collected from the Nangay and Jodrell 




Figure 1: Phase histogram of the Crab Pulsar obtained with photons 
above 100 MeV. The main radio peak is at phase 0. Two cycles are 
shown. 



radiotelescopes, that are involved in the LAT pulsar timing 
campaign The Crab ephemeris has been built with 
the TEMP02 fd\ package, using a set of 488 observations 
at 600 MHz and 210 observations at 1.4 GHz, collected 
between August 2008 and April 2009 [i,|3|. 
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Figure 2: Spectral energy distribution of the Crab Nebula in the energy range from 100 MeV to 300 GeV. The black continuous line is 
the fitted spectrum of equation [T] The dotted and dotted-dashed lines represent respectively the synchrotron and the Inverse Compton 
components of the spectrum. The spectral points have been obtained performing maximum likelihood fits in individual energy bins. The 
black error bars are statistical errors, while the grey error bars represent statistical errors added in quadrature to systematic errors. The stars 
are the EGRET results, and have been taken from ref. |12)I . 



3. Analysis of the Fermi LAT data 

3.1. Event .selection 

For this analysis 0, 01 we have considered a data 
sample covering the period from August 2, 2008 to April 
7, 2009. In particular, we have used only "diffuse" class [l[ 
photon events, that ensure the most effective background 
rejection. Events with zenith angles larger than 105° 
have been excluded from the analysis, in order to avoid 
contaminations from the Earth albedo gamma-rays. 

Figure [T] shows the Crab Pulsar phase histogram 
obtained selecting photons with energies above 100 MeV 
and with an angular separation from the nominal source 
position less than 

Oma. - max {[6.68° - 1.76° log^o E{MeV)] , 1.3°} 

where E is the photon energy. The phase corresponds 
to the position of the main radio peak. The Crab phase 
histogram exhibits two peaks at phases 0i = 0.9915 ± 
0.0005 and (/)2 = 0.3894 ± 0.0022 respectively The off- 
pulse window is defined in the phase range from 0.52 to 
0.87. 



3.2. Spectral analysis of the Nebula 

The spectral analysis has been performed following two 
different and complementary strategies. In the first case, 
the spectra are reconstructed using a maximum likelihood 
approach, that is implemented in the Fermi science tool 
''gtlike" [8|. In the second case, the spectra have been 
evaluated using an iterative unfolding procedure based on 
the Bayes theorem [13] . 

The maximum likelihood spectral analysis of the 
nebula has been performed selecting photons in the off- 
pulse window, in a 20° region centered on the nominal 
pulsar position, and with energies in the range from 
100 MeV to 300 GeV. The galactic background component 
has been modeled using GALPROP ll|, while the 
extragalactic component and the instrumental background 
have been described with a single isotropical power law 
spectrum. The Crab Nebula spectrum is well described 
by a sum of two power law spectra: 
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where Eq = 1 GeV. The two terms at the right hand side 
of eq. [1] are labeled with the indices S and IC because 
they are identified respectively with the falling edge of the 
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Figure 3: Spectral energy distribution of the Crab Pulsar in the energy range from 100 MeV to 300 GeV. The black continuous line is the 
fitted spectrum of equation [2] The spectral points have been obtained performing maximum likelihood fits in individual energy bins. The 
black error bars are statistical errors, while the grey error bars represent statistical errors added in quadrature to systematic errors. The stars 
are the EGRET results, and have been taken from ref. [I2l . 



synchrotron component and with the rising edge of the 
Inverse Compton component of the spectrum. 

The best fit values of the parameters are ks — (9.1 ± 
2.1 ± 0.7)cm-2s-iMeV-i, Ts = 3.99 ± 0.12 ± 0.08, 
kic = (6.4 ± 0.7 ± 0.1) cm-2s-iMeV-i and Tic = 1-64 ± 
0.05 ± 0.07, where both statistical and systematic errors 
are quoted. The latter, that are originated from the 
uncertainties in the LAT response, have been evaluated 
using modified instrument response functions with respect 
to the nominal ones. 

Figure [2] shows the spectral energy distribution 
E'^dN/dE of the Crab Nebula. The spectral points 
superimposed to best fit function have been obtained 
performing a maximum likelihood fit in individual energy 
bins. In this case a simple power law shape has been 
assumed for the source spectrum in each energy bin. In 
figure [5] the Fermi results are also_compared with the 
spectral point obtained by EGRET [1^. At energies below 
1 GeV the fiux measured by Fermi is higher than the 
one measured by EGRET. On the other hand, at higher 
energies, EGRET measured a higher fiux than Fermi . 

In the unfolding analysis, the Crab Nebula spectrum 
has been reconstructed without assuming any spectral 
model and without performing any fit. This analysis 
technique also allows to take into account the energy 



dispersion introduced by the detector response, that is 
not considered in the maximum likelihood approach. The 
observed count spectrum of the Nebula, used as input for 
the unfolding analysis, is built selecting the photons in the 
off-pulse window, and subtracting the background counts, 
evaluated using the background model obtained from the 
"gtlike" fit. The unfolded spectrum of the Nebula is in 
good agreement with the one obtained from the likelihood 
analysis. 

3.3. Spectral analysis of the Pulsar 

The spectral analysis of the Crab Pulsar has been 
performed with the same two techniques used for the 
Nebula. 

The maximum likelihood analysis of the pulsed 
emission has been done over the whole phase interval, 
including the Nebula in the background. The values of the 
spectral parameters of the Nebula reported in section 
have been renormalized to match the whole phase interval 
and have been fixed. The pulsar spectrum has been fitted 
with an exponential cutoff power law spectrum given by: 

dN / E\^^ / E\ 
where Eq — 1 GeV. 
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The best fit values of tlie parameters are k = (2.36 ± 
0.06 ± 0.15) cm-^s-iMeV-i, T = 1.97 ± 0.02 ± 0.06, 
and Ec = (5.8 ± 0.5 ± 1.2) GeV, where both statistical 
and systematic errors are quoted. We have also tried to 
perform a fit using a power law with a super-exponential 
cutoff, but the likelihood value was not significantly better 
than the one obtained using the function in eq. [2] 

Figure [3] shows the spectral energy distribution 
E'^dN/dE of the Crab Pulsar. The spectral points 
superimposed to best fit function have been obtained 
performing a maximum likelihood fit in individual energy 
bins, with the same procedure described in section 13.21 
In figure [3] the Fermi results are also compared with the 
spectral point obtained by EGRET 12]. The Fermi results 
are consistent with the EGRET ones in the energy range 
up to 8 GeV. 

The spectrum of the Crab Pulsar has been also 
reconstructed using the unfolding analysis. In this case the 
observed count spectrum of the Pulsar is built selecting 
photons in the on-pulse phase window, and subtracting 
the background counts in the off-pulse window, properly 
rescaled for the on/off phase ratio (i.e. 0.65/0.35 according 
to the definition in section 13. ip . As for the Nebula, the 
unfolded spectrum of the Pulsar is consistent with the one 
evaluated from the maximum likelihood analysis. 

4. Conclusions 

We have measured the energy spectra of the Crab 
Pulsar and Nebula in the energy interval from 100 MeV to 
300 GeV using a sample of data taken by the Fermi LAT 
during its first 8 months of operation in survey mode. 

The spectrum of the Nebula is well modeled with 
a sum of two power laws, with spectral indices = 
3.99 ± 0.12 ± 0.08 and F/c = 1.64 ± 0.05 ± 0.07, which 
describe respectively the falling edge of the synchrotron 
component and the rising edge of the Inverse Compton 
component. The pulsed emission is described by a power 
law spectrum with an exponential cutoff with spectral 
index F = 1.97 ± 0.02 ± 0.06 and cutoff energy E^ = 
(5.8 ±0.5 ±1.2) GeV. 

The spectral analysis has been performed using two 
different and complementary strategies, based respectively 
on a maximum likelihood fit and on an unfolding 
analysis. The spectra reconstructed with these two 
analysis techniques are in good agreement, thus confirming 
the reliability of our measurement. 

The Fermi LAT has allowed a precise measurement of 
the spectra of both the Crab Pulsar and Nebula that has 
been made possible by the wider energy range and by the 
better energy sensitivity with respect to its predecessors. 
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